?^U>-AOei  616 
UNCLASSIFIED 


NAVAL  POSTGRADUATE  SCHOOL  MONTEREY  CA  F/6  17/2*1 

CHANNELS  THAT  COOPERATIVELY  SERVICE  A  DATA  STREAM  AND  VDlCE  MES— ETCCU) 
NOV  79  DP  GAVER.  J  P  LEHOCZKY 
NPS55- 79-027 


NL 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  Of  THIS  PAGE  (When  Data  Enter e<t) 

f  REPORT  DOCUMENTATION  PAGE 


\  '■  report  auMSfea--, 
/  NPS55-79-jfc27 


2.  GOVT  ACCESSION  NO. 


14.  title  fend  ■ 


Channels  that  Cooperatively  Service  a  Data 
Stream  and  Voice  Messages*  I.  _ 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

3  RECIPIENT'S  CATALOG  NUMBER 


V  TYPE.OE  EtPWIWEWBB  By  VEREI 

)  Technical  ^ 

s.  PERFORMING  ORG.  REPORT  NUMBER 


it.  authorc*; 


D.  Py'fcaver  J.  P.^Lehoczky 


ft  CONTRACT  OR  GRANT  NUMBERr«J 

DD  1473 


[9  performing  organization  name  and  address 


Naval  Postgraduate  School 
Monterey,  Ca.  93940 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


.  '  61152N{RQO0g/il^l6 1 

'jjvfr  Q ,}.)  ;•  f  j  NO 00 14 80WR00 053 


11.  controlling  office  name  and  address 

Chief  of  Naval  Research 
Arlington,  VA  22217 


3.  NUMBER  OF  PAGES 


|  U  MONITORING  AGENCY  NAME  6  ADDRESSfl/  different  from  Controlling  Office)  I  15.  SECURITY  CLASS,  (of  thie  report) 


r/UL'~ 


Unclassified 


1S«.  DECLASSIFICATION,  DOWNGRADING 
SCHEDULE 


|  16.  DISTRIBUTION  STATEMENT  (of  this  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  it  different  from  Report) 


IS.  SUPPLEMENTARY  NOTES 


1 19.  KEY  WORDS  (Contfnue  on  reverse  aide  if  necessary  end  identify  by  block  number) 


Communication  Data-Voice  Communications 

Channels 

Queues 

Probability  Models 

20.  ABSTRACT  (Continue  on  reverse  aids  if  neceaaery  and  Identify  by  block  number) 

A  system  of  channels  mutually  accommodate  both  data  and  voice  messages, 
voice  having  pre-emptive  priority  but  being  a  loss  system,  and  data 
being  allowed  to  queue.  Approximations  to  the  data  queue  properties 
are  derived. 


DD  I  JAN  73  1473  EDITION  OF  1  NOV  65  |$  OBSOLETE 

S 'N  0 102-014- 6601 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS VaSF (Whan  Data  Entered) 

-3.54  yrn 


I .  INTRODUCTION 

A  system  of  channels  cooperatively  services  both 
voice  and  data  messages  arriving  at  one  node  of  a  communications 
network.  This  paper  is  devoted  to  the  analysis  of  a  particular 
channel-sharing  strategy,  in  which  voice  traffic  always  occupies 
its  channels  when  available,  but  data  service  is  allowed  to 
occur  on  empty  voice  channels.  Voice  traffic  is  taken  to 
be  of  high  priority;  voice  arrivals  that  find  all  voice  channels 
busy  are  treated  as  losses.  Note  that  voice  traffic  will 
be  relatively  infrequent  as  compared  to  data,  and  will  also 
exhibit  relatively  long  holding  (service)  times.  Data  traffic 
is  taken  to  be  heavy,  and  exhibits  very  short  holding  times 
(per  word  unit)  ;  compared  to  voice,  data  appears  to  arrive 
nearly  continuously;  when  all  data  (and  empty  voice)  channels 
are  filled,  queueing  occurs . 

k 

Research  in  part  sponsored  by  ONR  at  Naval  Postgraduate  School, 
N001480WR00067 ,  and  in  part  by  NSF  at  Carnegie-Mellon  University, 
ENG79  05526. 


1 


■  "**wi 


We  present  an  analysis  of  the  performance  of  a  special 
type  of  integrated  circuit  and  packet-switched  multiplexor 
structure.  This  structure  essentially  occurs  within  the 
SENET  network;  descriptions  are  given  by  Coviello  and  Vena 
(1975) ,  and  Barbacci  and  Oakley  (1976)  .  In  this  network  a 
time-slotted  frame  is  utilized;  a  certain  portion  of  each 
frame  is  allocated  to  voice  traffic,  while  any  remaining  data 
traffic  can  use  all  remaining  capacity,  including  that  left 
unused  by  voice.  Voice,  on  the  other  hand,  cannot  use 
capacity  unused  by  data,  and  operates  on  a  loss  system.  The 
subject  of  our  analysis  has  the  same  qualitative  flavor. 

Typical  performance  measures  that  may  be  calculated  are 
(i)  the  loss  rate  of  voice  traffic,  and  (ii)  the  expected 
waiting  time,  or,  equvialently ,  mean  queue  length,  of  the 
data . 

The  analysis  begins  with  standard  probabilitistic 
assumptions.  Specifically,  voice  traffic  arrives  according 
to  a  Poisson (A)  process,  and  each  customer  has  an  independent 
exponential ( n)  service  time.  Data  arrivals  are  according 
to  an  independent  Poisson (5)  process,  and  exhibit  exponential 
(n)  service  times.  A  total  of  c  channels  are  reserved 
for  exclusive  used  of  data,  while  v  channels  can  be  used 
by  both  data  and  voice;  however,  voice  pre-empts  data.  Voice 
operates  as  an  M/M/v  loss  system,  and  the  well-known  "Erlang  B" 
loss  formula  will  give  the  loss  rate.  We  are  mainly  interested 
in  the  behavior  of  the  data  queue;  however,  we  wish  to 
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develop  expressions  for  mean  queue  lengths  for  certain  extreme 

(and  realistic)  parameter  values.  First,  we  will  require 

that  6 /n  =  >  c.  This  assumption  indicates  that  the  data 

must  be  able  to  use  excess  voice  capacity  in  order  to  remain 

stable.  Second,  we  require  that  n/y  be  large,  perhaps  on 
4 

the  order  of  10  .  This  indicates  that  the  voice  requires  long 

service  periods  while  the  data  service  periods  are  very  short. 

This  problem  has  been  studied  in  a  number  of  papers 

including  Halfin  and  Segal  (1972)  ,  Halfin  (1972) ,  Fischer 

and  Harf.is  (1976)  ,  Bhat  and  Fischer  (1976),  Fischer  (1977), 

and  Chang  (1977) .  Many  of  these  studies  begin  with  the 

Kolmogorov  forward  equations  appropriate  for  this  system  and 

introduce  some  approximations  leading  to  a  solution.  While 

this  approach  is  entirely  appropriate,  the  approximations 

heretofore  introduced  have  not  been  tailored  to  the  p.  >  c 

d 

4 

and  n/y  ~  10  situation.  In  fact,  several  of  the  approxi¬ 
mations  give  quite  misleading  results  in  this  case.  We 
develop  an  approximation  which  is  tailored  to  these  rather 
extreme  but  realistic  parameter  conditions.  While  the 
Kolmogorov  equations  can  be  easily  written  for  the  Markov 
chain  {(Q(t),  N(t)),  t  >  0}  where  Q(t)  represents  the 
number  of  data  messages  in  the  system  and  N(t)  represents 
the  number  of  voice  messages  at  time  t,  the  fact  that  both 
Q ( t)  and  N(t)  are  subject  to  random  fluctuation  seems  to 
make  any  direct  approach  to  solving  the  equations  difficult. 

We  thus  propose  to  treat  the  data  as  a  deterministic  process 
behaving  like  a  fluid  flow. 


3 


2 .  THE  APPROXIMATION 

To  better  understand  the  behavior  of  the  data  queue 

process,  we  first  consider  the  special  case  c  =  0,  v  =  1. 

In  this  case,  the  data  can  use  the  single  channel  only  when 

voice  traffic  is  not  present.  Consider  a  set  of  parameter 

values  given  by  A  =  .01,  u  =  .01,  6  =25,  and  n  =  100.  It 

follows  that  p.  =  .25  and  p  =1.  The  overall  traffic 
d  v 

intensity  parameter  p  =  p^  +  p^d-q)  where  q,  the  voice 
blocking  probability,  is  .75,  so  the  system  is  stable.  Never¬ 
theless,  very  long  data  queues  will  occasionally  be  created 
for  the  following  simple  reason.  When  no  voice  traffic  is 
present  and  no  data  queue  is  present  the  system  appears  to 
data  traffic  to  be  an  M/M/1  system  with  p  =  .25.  There  will 
be  essentially  no  queueing  at  all,  and  this  situation  will 
persist  for  an  average  of  1/A  =  100  time  units.  However, 
when  a  voice  message  arrives,  the  channel  becomes  unavailable 
to  the  data,  and  all  data  messages  must  now  be  queued.  This 
queue  will  grow  at  a  rate  of  25  per  unit  time,  furthermore, 
the  voice  message  exhibits  a  long  holding  time  (on  the  average 
100  time  units) ,  so  the  data  queue  will  reach  a  height  of 
2500  on  the  average  before  it  can  begin  to  be  serviced.  The 
channel  is  now  free,  and  will  remain  so  for  about  100  time 
units — but  now  the  queue  has  2500  customers,  not  zero  as 
before.  It  is  clear  that  the  steady-state  mean  queue  length 
is  very  large  (2500  in  fact);  however,  it  is  also  clear  that 
this  classical  performance  evaluation  measure  can  be  very 
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misleading.  The  actual  behavior  of  the  data  queue  is  one  of 
long  periods  of  essential  emptiness  followed  by  long  periods 
of  great  queue  length.  The  mean  gives  an  average  of  these 
two  extreme  situations  and  therefore  is  misleading.  We  propose 
to  develop  approximations  for  this  mean  but  to  also  provide 
other  descriptions  of  system  behavior  such  as  idle  and  busy 
period  lengths,  first-passage  times  and  steady-date 
distributions . 

The  mean  queue  length  has  been  calculated  exactly  for 
c  =  0,  v  =  1  by  Fischer  (1977)  and  is  given  by 


(1+PV) 


(1-p) 


{II 

y 


pv+u  +  pvr} 


(2.1) 


where  p  =  +  pv/(1  +  pv).  It  is  clear  that  if  n/u  ~  104 

and  pv  ~  1  as  in  the  above  example,  then  of  the  two  terms 
in  brackets  n/u  will  be  large  compared  with  (1  +  p^)  ,  hence, 
we  can  ignore  this  term.  Ignoring  this  term  is  equivalent  to 
ignoring  the  queueing  that  occurs  when  the  system  is  empty. 

The  analysis  presented  in  this  paper  ignores  terms  of  this 
type. 

The  fluid  flow  approximation  is  based  on  treating  the 
data  as  a  deterministic  stream.  Let  us  suppose  that  there  are 
i  voice  channels  occupied.  This  leaves  c  +  v  -  i  available 
for  data.  Data  arrives  at  rate  6  and  is  serviced  at  rate 
(c  +  v  -  i)n  giving  an  overall  change  in  the  queue  length  of 
6  -  n (c  +  v  -  i)  =  r^  per  unit  time,  where  i  =  0,1,..., v.  It 
is  clear  that  rn  <  r.  <  • • •  <  r  .  We  assume  rrt  <  0  and 
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rv  >  0.  The  first  is  necessary  for  system  stability  while  the 
latter  follows  from  >  c.  Thus  there  is  a  state  N  for  which 

rN  £  0  C  rN+1.  Wt  treat  the  case  rN  <  0 ,  while  rN  =  0  is 
a  straightforward  generalization.  We  refer  to  the  states 
0,  1,  ...  ,  N  as  "down"  states,  while  N+l,  .  ..,  v  are  "up" 
states  (0  <_  N  <_  v,  so  the  two  sets  of  states  are  nonempty)  . 

These  names  reflect  the  fact  that  if  i  voice  channels  are 
occupied,  then  the  data  queue  tends  to  increase  if  i  is  an 
up  state,  and  to  decrease  if  i  is  a  down  state.  The 
steady-state  distribution  of  the  occupancy  of  the  voice 
channels  is  given  by  a  truncated  Poisson (p^)  distribution, 


0  <  i  <  v 


(2.2) 


and  the  loss  probability  q  =  p  . 

I 

For  the  data  queue  to  remain  stable 


I  r.pj/i!  <  0  . 
i=0  1  v 


If  one  defines  p  =  +  p  (l-q)]/(v  +  c)  then  the  stability 

condition  becomes  p  <  1. 

We  wish  to  compute  a  variety  of  quantities  for  the  data 
queueing  system.  These  quantities  include 


Hi (x) 


P (voice  is  in  state  j  when  queue  empties  I  voice 
is  in  state  i  and  data  in  state  x) ,  0  <  i  <  N, 
0  <  i  <  v.  _  - 
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x^(x)  =  expected  first-passage  time  for  data  from  state  x 
to  state  0  starting  in  voice  state  i,  0  <_  i  <_  v. 

a.(x)  =  expected  area  under  data  queue-length  process 

accumulated  during  the  first  passage  time  to  0, 

0  _<  i  <_  v. 

The  above  quantities  give  important  characterizations 

of  the  system  performance.  The  first-passage  times  indicate 

the  time  needed  to  work  off  a  backlog  of  size  x.  The  area 

gives  essentially  the  waiting  time.  If  the  queue  is  empty  and 

the  voice  is  in  a  down  state  then  for  the  fluid  model  the  queue 

will  remain  empty  until  the  voice  reaches  the  first  up  state, 

N+l.  The  queue  immediately  begins  to  grow  at  rate  rN+i‘  It 

follows  that  tn+i ^  represents  the  expected  duration  of  the 

busy  period.  Similarly  aN+i^  gives  the  expected  area 

accumulated  during  the  busy  period.  Using  renewal-theoretic 

ideas  aN+-^  (0) /tn+^  (0)  gives  the  mean  queue  length  during  the 

busy  period.  Similarly  pN+^  j_(0)  gives  the  probability  that 

the  busy  period  will  end  in  voice  state  i.  The  time  for  the 

voice  to  reach  N+l  from  i  and  hence  the  expected  time  to 

initiate  a  new  busy  period  is  easily  calculated  from  the 

birth-death  nrocess.  Let  us  designate  this  mean  by  S..  Then 
rN 

ii_Q  S^Pn+i  j_(0)  =  T  gives  the  expected  idle  time  (we  ignore 
all  queueing  during  this  period).  Clearly  T/(T  +  TN+^(0)) 
gives  the  steady  state  data  component  idleness  probability  and 


m 


aN+l^)"T/,^T  +  tN+1  •  TN+^(0)  gives  the  steady-state  mean 

data  queue  length.  It  is  clear  that  the  quantities  a  (x)  , 
xi(x)/  and  j ( x)  give  valuable  insight  into  the  behavior 
of  the  queueing  process,  incidentally  providing  all  of  the 
standard  queueing  performance  measures.  The  special  case 
of  one  down  state  (N  =  0)  is  easiest  to  handle,  since  in 
this  case  the  P^j's  can  be  ignored. 
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3.  DERIVATION  OF  p_(x)  FUNCTIONS 

We  use  a  backward  equation  approach.  Let  us  assume  that  at 
time  t  =  0  the  queue  length  is  x  >  0  and  i  voice  channels 
are  occupied.  It  follows  that  at  time  dt,  the  new  queue  length 
will  be  x  +  r^  dt.  The  system  will  remain  in  state  i  with 
probability  1  -  (X  min(l,v-i)  +  iy)dt  +  o(dt),  will  move  to 
state  i+1  with  probability  X  min(l,v-i)dt  +  o(dt),  or  will 
move  to  i-1  with  probability  iydt  +  o(dt).  Thus 

p^ j  (x)  =  p^_,  (x  +  r^dt) (1  -  (X  min(l,v-i)  +  iy)dt  +  o(dt)) 

+  p.  .  . (x  +  r.dt)  iy  dt 

l-l  j  l 

+  P^+1  j (x  +  r^dt) X  min(l,v-i)dt  +  o(dt)  (3.1) 

One  expands  the  (x  +  r^  dt)  into  p^j  (x)  +  r^p.V  (x)dt 

+  o(dt),  collects  terms  and  lets  dt  0  to  derive 

p! . (x) (-r.) 

=  -p .  . (x)  ( X  min ( v-i , 1 )  +  iy)  +  p.  .  (x) iy 

1 J  1  -*■  r  D 

+  P^+1  j  (x)  Xmin(v-i,l)  with  0  <_  j  N.  (3.2) 

If  r^  =  0 ,  then  (3.2)  indicates  a  linear  relationship  among 
PNj(x),  PN-1  j(x)  and  PN+1  j(x).  This  relationship  serves  to 
allow  elimination  of  p^_.(x)  and  therefore  allows  us  to  assume 
r^  <  0 .  Equation  (3.2)  can  be  divided  by  -r^  and  the  entire 
system  rewritten  in  matrix  form  to  yield 
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P'  (x)  =  ^  Q*  P(x) 


here  P(x)  =  (p„  (x)  )  ,  0  <_  i  ^  v,  0<j£N,  a  (v+l)x(N+l) 

* 

stochastic  matrix  for  each  x,  and  Q  is  defined  by 


where  k  =  c  +  v  -  p,.  We  assume  0  <  k  <  v. 

d 

Equation  (3.3)  can  be  routinely  solved  to  aive 

P(x)  =  exp(^  Q* )  P(0)  (3.4) 

2 

where  exp(M)  =  I  +  M  +  M  / 2 1  +  •••  for  a  square  matrix  M. 
Interestingly,  one  still  needs  to  determine  P(0)  before 
P(x)  is  fully  determined. 

To  determine  P(0),  we  partition  into  down  and  up  states. 

Thus 

/  Vx)  \ 

P(x)  = 

PUU)  / 


where 
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3.  DERIVATION  OF  p..(x)  FUNCTIONS 

We  use  a  backward  equation  approach.  Let  us  assume  that  at 
time  t  =  0  the  queue  length  is  x  >  0  and  i  voice  channels 
are  occupied.  It  follows  that  at  time  dt,  the  new  queue  length 
will  be  x  +  r^  dt.  The  system  will  remain  in  state  i  with 
probability  1  -  (X  min(l,v-i)  +  ip)dt  +  o(dt),  will  move  to 
state  i+1  with  probability  X  min(l,v-i)dt  +  o(dt),  or  will 
move  to  i-1  with  probability  ipdt  +  o(dt).  Thus 

p^ j (x)  =  p^j (x  +  r^dt) (1  -  (X  min(l,v-i)  +  ip)dt  +  o(dt)) 

+  p.  .  . (x  +  r.dt)  ip  dt 

*1-1  j  i 

+  j (x  +  r^dt) X  min(l,v-i)dt  +  o(dt)  (3.1) 

One  expands  the  p-.(x  +  r.  dt)  into  pi  •  (x)  +  r.p.'.(x)dt 
+  o(dt),  collects  terms  and  lets  dt  -*•  0  to  derive 

P| j (x) (-r.) 

=  -p .  . ( x)  ( X  min(v-i,l)  +  ip)  +  p.  ,  .  (x) ip 

+  Pj_+i  j  (x)  Xmin(v-i,l)  with  0  £  j  <  N.  (3.2) 

If  rN  =  0 ,  then  (3.2)  indicates  a  linear  relationship  among 
PNj(x),  pN_1^ j (x)  and  PN+1  j (x) .  This  relationship  serves  to 
allow  elimination  of  pN_.(x)  and  therefore  allows  us  to  assume 
r^  <  0 .  Equation  (3.2)  can  be  divided  by  -r^  and  the  entire 
system  rewritten  in  matrix  form  to  yield 
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PD(x)  =  (p±j  (x)  )  ,  0  <  i  <_  N,  0  <  j  <  N 

P  (x)  =  (p..(x)),  N+l  <i<v,  0  <  j  <  N  . 

~U  l]  —  —  —  — 

If  no  data  queue  is  present  (x  =  0)  and  i  is  a  down 
state,  then 

f  +  1  if  i  =  j 


since  emptiness  is  instantaneously  achieved-  Thus  p^OD)  =  I» 
the  (N+l)  x  (N+l)  identity  matrix.  It  remains  to  calculate 

V01- 

A  second  system  of  equations  can  be  developed  as  follows. 
Beginning  in  state  i  (i  up)  and  x  >  0,  one  must  first  return 
to  level  x  and  then  hit  0.  The  return  to  level  x  must  occur 
in  a  down  state.  This  allows  one  to  write  a  system  of  "Chapman- 
Kolmogorov  like"  equations 


P#U(X)  =  Py  ( 0  )  PD(X) 


(3.5) 


Equations  (3.3)  and  (3.5)  can  be  combined  to  aive  an  exoression 
for  P^fO).  This  expression  is  in  the  form  of  a  matrix  quadratic 


equation : 


<Pq(0> ' 


N+l 


pu(0) 


0 


(3.6) 
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This  equation  can  be  rewritten  to  be 


Pj.  (0 )  A.  . 

— u  ^1 1 


^21 


~u  ^ 


^22?U(0)  =  & 


(3.7) 


where 


with 

A-^  an 

(N+l) x (N+l ) 

matrix 

and 

&22  a 

( v-N) x ( v-N) 

matrix 

Equation  (3.7)  does  not  yield  a  closed  form  solution  except 
in  very  special  cases.  It  can,  however,  be  solved  numerically 
using  a  Newton-type  iteration.  Such  solutions  have  been  carried 
out,  but  the  results  will  not  be  provided  here. 


DERIVATION  OF  (x)  FUNCTION 

The  first-passage  time  functions  can  also  be  derived  using  a 
backward  equation  approach.  A  straightforward  derivation  gives 
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where  p^OD)  has  been  previously  determined.  Using  the  partitioned 
version  of 
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Equation  (11)  becomes 


V0) 


y  (~U(0)~12 


i22)'1  <-Pu(0),-Iv.N) 


/  r  \ 


u/ 


(3.12) 


For  example,  in  the  special  case  v  =  1,  PyfO) 


=  (1), 


k  =  1-Pj  and 


xu(o)  x ! ( o )  -  a  -  p)  • 


(3.13) 


In  general  if  N  =  v-1  (only  1  up  state) 

v-1 

l  p  .  (0)  - i — 

1  i=0  VI  K-l  K-V 

T  (0)  =  -  i— _ 

v  y  p 


K-(V-l)  PV, V-l(0^  +  K-V 


(3.14) 


while  if  N  =  0  (only  1  down  state) 


rn(0)  =  -  -±- 
~U  Ky 


f1' 

\\u 


i-1 


(Pv  0-0)  -  a22 


/A 


\ 


k-2 


K-V 


(3.15) 


Explicit  solutions  for  other  cases  can  be  written  down  but 
become  very  complicated. 
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4.  DERIVATION  OF  a^x)  FUNCTIONS 

The  backward  equation  approach  gives  a  straightforward 
derivation  of  the  area  accumulated  under  the  queue  length  process 
during  the  first-passage  time.  One  derives 


Clearly  aD(0)  =0,  but  a^CO)  must  be  determined.  Equation 
(4. 1)  can  be  rewritten  as 


which  has  a  straightforward  exponential  solution  once  the  initial 
conditions  have  been  determined.  To  this  end,  a  second  set  of 
equations  can  be  found  using  the  Chapman-Kolmogorov  approach. 


Beginning  in  an  up  state  at  level  x  the  process  must  return 
to  level  x,  then  to  0.  The  expected  area  accumulated 
during  the  return  to  x  is  given  by  a^(Q)  +  xt^(0).  It 
follows  that 


a  (x)  =  au(0)  +  x  t  u  ( 0 )  +^(0)  a(x).  (4.3) 


Equations  (4.2)  and  (4.3)  can  be  combined  to  give 


-’b101  -  S  (pul0) •'Iv-N>  0 


3y  partitioning  Q  one  finds 


V0) 


(4.4) 


V°>  = 


-  S’  'V0)%2  ~  $22>  TU(0) 


(4.5) 


A  simple  example  of  the  calculation  involved  in  (1.5) 
is  the  v  =  1  case .  It  can  be  shown  that 


p,(l-pj 
/n  \  n  ^  ^ 

a  (0)  =  —y  - 2 - 

1  u2  (1  +  Pv)2(l-P) 


for  p  <  1 


(4.6) 


Carrying  this  a  sten  further  one  can  see  that  ths 
idle  period  has  mean  length  1/A  =  (1/u) (1/p^) .  Thus 
recalling  (3.13),  we  find  the  mean  queue  length  given  by 
a1(0)/(r1(0)  +  1/A)  or 


e(q)  =  a 


p  <  1  . 


(4.7) 


pvpd 


(1-P) (1+PV) 


The  v  =  1  case  can  be  carried  further.  Once  a,(0) 
and  1,(0)  are  known,  a,  (x)  and  i^fx)  can  be  determined. 
Equation  (3.11)  becomes  t^x)  =  t1  (0 )  +  tq(x)  indicating 
x,(x)  -  tq(x)  =  t,(0)  given  by  (3.13) .  Equation  (3.9)  can  be 
routinely  solved  to  find 


ro(x)  "  n(3Tp^r  (1  +  (1+pv)  (i-p)  )x 


(4.8) 


Ti(x)  =  - 


U 


H  V  l 

nTI-^^T  (1  +  (1+pv)  (l-p))x  +  {(1+pv)  (1-p)' 


The  area  function  can  also  be  explicitly  determined. 
Equation  (4.3)  becomes  a,(x)  -  a^(x)  =  xt,(0)  +  a,(0). 
Substituting  this  into  (4.1)  gives 

a„(x)  -  %-  (i  +  ipO))  +  £  S.  ai(0)x 

(4.9) 

al(x)  =  \  <7  +  T1(0))  +  £  3,(0)  +  r1(0))  x  +  3,(0), 

all  coefficients  of  which  have  been  previously  determined. 

In  the  special  case  =  1/4,  pv  =  1,  p  =  3/4  mentioned 

earlier  with  A  =  .01,  u  -  .01,  6  =  25,  n  =  100  then 
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T0(xl 

(x)  =  4x  +  200 
a^(x)  =  x2  +  lOOx 
ax(x)  =  ^  x2  +  300x  +  750,000 


E(Q) 


a1(0) 

T.  (0)  +  -  — 

1  W  PV 


750,000 
200  +  100 


(4.10) 


2500  . 


The  voice  loss  rate  is  given  by  Ap  /(1+p  )  =  A/2 


.005. 


5.  STEADY-STATE  DISTRIBUTION  OF  DATA  QUEUE  LENGTH 

One  can  use  a  forward  equation  approach  to  develop  an 
equilibrium  distribution  for  the  data  queue  length.  Define 
p(x,j,t)  to  be  the  probability  of  j  voice  channels  occupied 
and  x  data  units  in  the  system  at  time  t.  It  is  easily 
seen  that  for  x  >  0  and  dt  small 

p(x,j,t  +  dt) 

=  p  (x-r  jdt ,  j  ,  t)  (1  -  (A  +  jy)dt)  +  pfx-r^^dt,  j-l,t)Adt 
+  p  (x-r_.+1dt,  j+1 ,  t)  (j+1)  y  dt  +  o  ( dt)  ,  0  <_  j  <  v 

where  p(x,-l,t)  =  0. 

Standard  manipulations  that  treat  x  as  a  continuous 
variable  lead  to  these  equations 

r.  +  =  -(A  min  (1 ,  v- j)  +  jy)  p(x,j,t) 

3  3t  (5.1) 

+  A p (x , j -1 , t)  +  (j  +  1) y  p(x,j+l,t),  x  >  0. 

Setting  t  -*■  »  and  assuming  p(x,j,t)  p^  (x)  and 
[  3p  (x,  j  ,  t)  ] /3 1  -*■  0,  we  find 

Tj  p  ^  (x)  =  -(A  min  (1 ,  v- j  )  +  jy)  p^.  (x) 

+  A  p^^tx)  +  ( j  +  1 )  y  p  (x).  0  <_  j  <  v  (5.2) 

with  x  >  0,  p  ^(x)  =  pv+^  (x)  =  0.  The  prime  denotes  x-derivatives . 
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Equation  (5.2)  is  incomplete  as  it  does  not  contain 
information  about  the  boundary  behavior. 

Equation  (5.2)  can  be  summarized  in  matrix  form  by 

P' (x)  =  -P  (x)  R*  ^  ,  x  >  0  (5.3) 

where  P(x)  =  (pQ(x),  ...  ,  p  (x))  and 


Equation  (3)  can  be  routinely  solved  to  get 

P(x)  =  c  exp(-  ^  R  x)  (5.4) 

with  c  =  (Cq ,c^ , . . . , c  ) .  The  constants  c  must  be  determined, 
and  equation  (5.4)  gives  only  the  density  function,  not  the  mass 
at  the  boundary.  In  view  of  the  fluid  flow  approximation, 
there  will  be  mass  at  0,  given  by  for  each  down  state  i, 

0  £  i  <_  N,  however,  no  mass  at  the  boundary  for  any  up  state, 

N+l  <.  i  <.  v.  Furthermore,  the  equilibrium  distribution  over  i 
is  given  by  ( c */i ! ) j 1  It  follows  that 
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This  gives  v+1  independent  equations  which  determine  c. 
Once  c  has  been  determined,  p^(x),  0  <_  i  £  N  are  determined 
by  (5.4).  One  can  now  determine  ,  . .  . ,  ttn  ,  the  boundary 
probabilities,  using  (5.5).  The  equilibrium  distribution  is  now 
completely  determined. 

Let  us  consider  the  special  case  v  =  1.  Equation  (1) 
becomes 

p- (x)  =  _  H  ^  po(x)  +  i  Pl(x)) 

Pl(x)  =  '  n  po(x)  '  pl(x))  (5-8) 

kp^fx)  +  (k-l)p|(x)  =  0, 
hence, 

kpQ(x)  +  (k-l)p1(x)  =  0  and  p-^x)  =  -  pQ(x). 

Substitution  into  (8)  yields 
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*0(x)  = 


y  ,"pv  l 


-  £  (  <  ‘  po(x)  = 


E  (l-P)  <i+pv> 

n  pd(1"pd) 


Pft(x) 


p0(x)  =  C  exp  [  J 


(1-p) (l+p„) 


Pd(1_pd] 


(5.9) 


p^x)  =  c 


(1-p)  (1+Pj 


d  Kd 


1+V  x] 

rp^r  XJ 


OO  p 

/  p1(x)ax  =  J-f2- 


thus 


(1"pd)  n  pd(1~pd} 


and  c  = 


pv(1~p)  y 


(1-Pd)  n 


(5.10) 


It  is  then  possible  to  determine  ttq  by  means  of  the 


relation 


oo 

"o =  i+p^:  _  l  po(x)dx 


Finally  we  find  the  marginal  queue  length  distribution  in  a  simple 


explicit  form 


jj 

n  <1'l’d)20d 


n  pd(1"pa) 


x  >  0 


p(x)  = 


(5.11) 


x  =  0 


The  mean  queue  length  is  given  by 


E(Q)  =  (£) 


rdpv 


(1-p) (1+PV) 


P  <  1.  (5.12) 
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6.  GENERALIZATIONS 

The  preceding  analysis  giving  the  data  queue  length 
can  be  carried  out  in  greater  generality.  One  may  wish  to  view 
voice  transmissions  as  a  stream  of  alternating  bursts  and  dead 
times.  At  any  time,  t,  a  voice  customer  assigned  to  a  particular 
channel  will  be  in  one  of  two  states:  active  transmission  or 
inactive.  During  an  inactive  period  the  channel  could  be  used 
for  data  transmissions.  Such  a  strategy  could  greatly  increase 
the  data  capacity,  or  reduce  the  data  queue  length,  or  both. 

One  might  assume  a  voice  user  moves  between  the  active  and  inactive 
states  in  a  Markovian  way.  As  a  result,  one  can  define  voice 
states  (1,  ...  ,  M>  and  a  continuous  time  Markov  chain  with 
generator  g  describing  the  movement  among  these  states.  For 
each  voice  state  i  the  number  of  channels  available  for  data 
transmission  can  be  found,  and  a  rate  of  increase  or  decrease  in 
the  data  queue  length  determined.  Let  that  rate  be  denoted  by 
and  assume  that  the  voice  states  are  labelled  in  such  a  way 
that  r-^  <_  r2  <_  •  •  •  <_  rM-  We  would  assume  there  is  a  state  I 
for  which  r^.  <  0  <  rJ+1  Up  and  down  states  can  now  be  defined. 
The  analysis  carried  out  in  previous  sections  can  be  shown  to 
hold  for  this  more  general  situation.  The  expression  ug*  must 
be  replaced  by  g,  while  the  expression 
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must  be  replaced  by 


It  follows  that  first  passage  times,  queue  lengths, 
and  busy  period  lengths  can  be  determined  for  this  more  general 
problem. 

A  second  generalization  involves  a  voice  limitation 
procedure.  One  way  to  prevent  the  buildup  of  very  long  data 
queues  is  to  increase  the  number  of  voice  channels.  This  can 
be  accomplished  by  providing  lower  transmission  rates  on  certain 
voice  channels,  for  example,  one  might  provide  10  8kBPS  channels; 
however,  when  8  of  those  are  in  use  one  might  divide  the  remaining 
2  into  4  2kBPS  channels.  This  results  in  lower  quality  trans¬ 
mission  but  less  delay.  Again  this  situation  can  be  modelled  by 
assuming  various  voice  states  {1,  ...  ,  N)  and  movement  among 
them  according  to  a  Markov  chain  with  generator  Q.  Each  voice 
state  determines  a  data  queue  rate  r^.  Again  the  preceding 
analysis  can  be  applied  directly  to  this  case.  One  can  there¬ 
fore  study  the  tradeoffs  between  data  queue  length  (and  delays) , 
voice  blocking  probabilities,  and  voice  transmission  quality 
from  formulas  in  this  paper. 
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